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SUMMARY 



This report gives a simple method of estimating the 
torsional stiffness of thin shells, snoh as- "box teams "or 
stressed- skin wings under large torque loads. A general 
efficiency chart for shells in torsion is established, 
based on the assumption that the efficiency of the web 
sheet in resisting deformation decreases linearly with 
the average stress. The chart is used to calculate the 
torsional deflections of eight hoz "beams, a test wing pan- 
el, and a complete wing; the results of the calculations 
are shown in comparison with the test results. The agree- 
ment is probably as good as might be expected considering 
the empirical nature of the method and the well-known dis- 
persion between results of tests on thin sheet-metal 
structures . 

INTRODUCTION 



A knowledge of tho torsional stiffness of box "beams 
and Btressed-skin wings under large torsional moments 1b~ 
necessary for calculating the interaction between sparB 
(refer once' 1).. Thin shells of this nature bucklo under 
small loads and tho covor platos work no longer in shear, 
but in diagonal tension. This transition decreases the 
torsional stiffness considerably, hut it 1b not sufficient 
to explain the large deformations that tests show, partic- 
ularly under large loads. The experimental lqad-aef orma^- 
tion diagrams are cuVes rather than straight lines, indi- 
cating that the torsional stiffn.es b decreases continuous- 
ly as the load increases. A nethod of correcting the 
theoretical- stiff ness constants was therefore developed 
based on the assumption that the efficiency of the sheet 
metal in resisting deformation decreases linearly with 
the average stresB in the sheet, the rate of decrease de- 
pending on the aspect ratio of the individual sheets. 
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CALCULATION 07 THE TORSIONAL STIFFNESS 07 SHELLS 



When concentrated torques T act on the ends of a 
tube (fig. l) of length L the resulting angle of" twist 
is 

• - & &> 



where 0 is the modulus of shear and 



J is the torsion constant of the section. 



For a thin-walled tube the torsion constant is 

J = 2±- .._ _( 2 ) 



where A is tho area bounded by the median line of the 
cross section, ds is a differential element of the per- 
imeter and t is the thickness of this element. The 
shearing stress at any point is 



Formulas (2) and (3) are derived under the assump- 
tions that the wall thicknesses are small and that the 
cross sections do not change their shape. Both formulas 
have been verified experimentally by various investigators 
and have been found to check the experimental values 
within very close limits for round, rectangular, and 
streamline tubing of the thickness-diameter ratios UBual 
in aircraft construction. 



Box spars and stressed- skin wings are also thin- 
walled tubes, but tho wall thicknesses are very much 
smaller in relation to tho dimensions of tho cross section 
than in tho case of tubing. Whon such a thin shell is 
subjected to torsion, the sides will buckle at a very 
small load and will be transformed into diagonal-tension 
field beams. 
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lor the present purpose it is convenient to calculate 
the deflections oocurring after the diagonal- tension field, 
has formed by using the formulas for shear deflection. 
Consider a square sheet with a side length of unity and a 
thickness t, subjected to a shearing force* S. If . -the 
sheet did not buckle, the shear dofloction would.be 




where •& is the modulus of shear. If the same sheet were 

a panel of a Wagner beam, i.e.', Tf IT did buckle,' the de- 
flection would be 

' - tf - <« 



E •'being the modulus of elasticity. When (4) and (5) are 
compared it will be seen that (4) can be UBed to calculate 
the deflection of a Wagner beam panel provided that the 
actual thickness t is replaced by an effective thickness 



or, for duralumin t e = 5/8 t, 

• formula (5) is based on the elementary theory of the 
Wagner beam presuming rigid flanges. In the particular 
type of structures under consideration, however, the 
flanges of the beams are very flexible. In some ca¥es the 
flanges consist of angles of a lesser, thickness than the 
web itself. lor beams with such flexible flanges, the ap- 
plicability of the theory 1b very doubtful, it is only 
known that the stress distribution in the web loses uni- 
formity more and more while the load increases, the stress 
concentrating along the folds as indicated in figure 2. 
The effect of this change in stress distribution on the 
shear deflection can be indicated by writing .. — """ 




where T| is an efficiency factor depending on the average 
stress and on the flexibility of the flanges. Since r\ 
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decreases with increasing' BtroBSv the tor_sional stiffneBS 
o^f s/.&hel'l is not a. constant , but. decreases with increas- 
ing torque , resulting". J.n a curved, load-deformation diagram. 

A convenient means of establishing empirical rela- 
tions for T) was afforded by the test results given in 
reference 3-', ' which describes aseries of tests on box beams 
of similar construction. The beams consisted of smooth 
sheets, 0,010 to 0.041 inch thiok, riveted to 4-inch chan- 
nels, 0.015- to 0.049 inch thick. (See figs. 5 to 12.) 
The bulkheads were made of sheet 0.037 to 0.049 inch thick. 

it was. found that T| could be assumed to decrease 
"linearly with increase of stress from unity at zero stress 
(fig. 3). The average shearing stress f fi calculated by 
formula (3) was used as an index for convenience, although 
this stress does not actually exist as a shearing stress 
after the diagonal-tension field has formed. The rate of 
Change of efficiency with stress waB found to be only a 
function of the spacing of the uprights supporting the 
flanges of the Wagner beams or* in other words, of the as- 
pect ratio d/h ' of the component plates (fig. 2), This 
fact was eomewhat surprising, because one would naturally 
expect the cross section of the flanges to have some in- 
fluence. For the very flexible flanges, however, which 
are a characteristic of the beams investigated, the evi- 
denoo appears to indicate that thi b influenoe is of a mi- 
nor nature. Undoubtedly, flanges consisting of heavy an- 
gles, would have some such influenoe. The beam of figure 
12 has corner reinforcements, but as the teBt was not car- 
ried t,o failure, no conclusions could be based on it. 

.. The procedure of calculating the torsional stiffnesB 
for a given beam or shell under any given torque will 
therefore be as follows: 

(l) Calculate the critical buckling stresses in 
shear for the component sheets. The edges of the- sheets 
should be assumed to be simply supported, because the 
beneficial effect of the existing elastic restraint is 
probably moro than canceled by tho dotrimontal offoct of 
initial buckles. Tho critical stross calculated is prob- 
ably not conservative but is, in general, sufficiently 
accurate for the present purpose. Figure 4 gives a graph 
for calculating the buokling stresses of duralumin sheets 
based on values given in reference 3. 
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(2) Calculate f 0 for each sheet by formula (3). 
If this stress 1b higher than the critical buokling stress 
found in step (l) , use the efficiency chart (fig. 3) to 
find T|. 



(3) Calculate the torsion constant . J by formula 
(2) , substituting for the thickness the actual value t 
if the sheet has not buckled, or the value t e X T| = 5/8 tf) 
if the sheet has buckled. 

In structurally complex shells, e.g., wings with 
several spars and longitudinal stringers, it may be dif- 
ficult to decide the proper value of h to use. If longi- 
tudinal seams exist between spars in the skin, as on the 
panel shown In figure 13, the distance between seams should 
be used. In any complex, built-up shell, values of d/h 
less than 0.8 should be used with caution and substantiated 
bytests if necessary. 

If it is desired to construct a load-deformation dia- 
gram, the load is divided into a number of .convenient in- 
tervals, say 2,000 in.rlh. She stiffness is then calcu- 
lated for the middle of each interval (1,000 in. -lb., 
3,000 in. -lb,, etc.) and the increment in angle of twist 
due to each increment of torque is calculated. The Incre- 
ments of twist are plottod against torque; the resultant 
broken line gives a series of tangents to the desired 
curve. 



COMPARISONS BETWEEN CALCULATED AND EXPERIMENTAL 
■ ■ DEEORMATION DIAGRAMS 



Eigures5 to 12 show the calculated deflection ..curves, 
and the experimental points for the beams of reference 2. 
The construction of all the beams is very 3imilar. The es- 
sential dimensions are given on each figure. j. 

The beams of figures 5, 6, and 7 are idontical with 
the excoptlon of bulkhead spacing. The boam of figure 9 
is the only one of the series 'that shows considerably more 
deflection than predicted.. No explanation could be found 
for this behavior. Comparative weight estimates indicate, 
however, that some sheets of this beam may have had less 
than nominal thickness. It must be remembered that the 
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commercial tolerance for thiokness variation is 0.0016 
-inch or 7 peroent of the top and bottom cover. 

The "beam of figure 10 has a cambered top cover. She 
buckling stress for this curved sheet was assumed to be 
the same as if the sheet were part of a cylinder and was 
calculated by the method given in reference 4. . . 

The results on the wing panel (fig. 13) were taken 
from reference 5. The sheets are supported by ribs and 
stringers bo that they are. divided into nearly square 
component panels. The effect of ..the curvature on the 
critical buckling stress was neglected;, i.e., it was as- 
sumed that the leading edge buckled at the same. time as 
the flat part of the covering; so that t Q = 5/8 t over 
the whole surface. Tai s assumption is contrary to the 
experimental evidence; tests always show that the sharply 
curved nose covering does not buckle until large loads 
are applied; in fact, buckling of the leading edge is 
usually the causo of ultimate failure of tha wing shBll 
in torsion. (See fig. 14..) However, the part of the 
nose that does not buckle is. a small part of the whole 
cover. Furthermore, the cover does buckle a short dis- 
tance behind the nose, where the curvature is still very 
largo; tho bucklos occurring at .this location tend to 
pull tho shoot down to chord linos subtending -the actual 
curve of the profile, thus reduolng- the area included by 
the profile forward of the front spar. Thi s reduction 
of included area reduces the stiffness (formula (2)); 
therefore the assumption of all the covering buckling 
tends to give an average. The fact that the actual stiff- 
ness was larger than the calculated- value is partly due 
to the fact that tho web of the front spar was noglootod 
in the calculation of the torsional stiffness. Inner 
walls, such as this spar, usually contribute less than 5 
percent to the torsional stiffness of a shell but in this 
case the more exact formula, which includes the effect of 
the front spar, gave a torsional stiffness 7 percent high- 
er than that calculated .by tho simple- formula (2). 

Figuro 14 shows tho calculated and tho oxporlmontal 
twist at 5 stations of an all-metal wing (reference 6). 
Ho data were available on the spar system except photo- 
graphs, which Bhowed a multispar system similar to the 
Junkers type (transposed flanges) with apparently very 
light web members. Calculations given in reference 1 
show that for 2-spar* stressed-skin wings the influence 
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of the spars on the torsional stiffness is negligible; 
the. same thing is probably true for this multispar wing. 
The value d/h = 1 was used, because it is representa- • 
tive of the major part of the cover. 

■ * 

Failure occurred between stations 4 -and 6 by buck- 
ling and rupture of the leading odge. The result was 
the large deformation recorded at. station 5 for the last 
two load increments. 

Omitting curve 2 of figure 14, because there was 
obviously some disturbance of the apparatus before the 
last three readings wore taken, and also curve 5 of tho 
same figure because failure occurred here, the average 
ratio of calculated to observed deflection for the last 
points of all test curves given 1b 0,92. 



CONCLUSION 



The method given for estimating the torsional stiff- 
ness of a thin shell under large . torques cannot lay claim 
to great accuracy* It is empirical in nature and based 
on not very extensive evidence. Additional experimental 
or theoretical research may replace tho straight linos of 
the efficiency chart by curves, change the spacing of the 
curves, and introduce the stiffnesses of. the flanges and 
of the bulkhead as factors. Until such additional work . 
is done, the chart may serve as a guide. 



Langley Memorial Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Held, Ya. , June 11, 1934. 
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Figure 1. 
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Figure 3.- Effioienoy chart for thin shells in torsion. 
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Figure 4.- Critical shearing BtreBB for duralumin. 
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Figure 5.- Computed and experimental Talues of variation of angle of twiit witt 
torque for box oeaae. Ixperimantal data from referenoe 8. 
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Figure 6.- Oomputed and experimental values of variation of angle of twist 
vita torque for 'beams. Experimental data from referenoe S. 
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Figure 7.- Oomputed and experimental values of variation of angle of twlit 
with torque for box beams, experimental data from referenoe 8. 
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Figure 8.- Computed and experimental values of variation of angle of twist 
with torque for box beams. Experimental data from referenoe 8. 



10,000 



8,000 



^ 6,000 



g 4,000 

E-i 



2,000 





















1 

Failure 


















* — 




X 
























X 


— X 






















f 


X 


X 
























- X 


< 


























X 


























X 

X 














n.naT 


t 










I x 










0 


.041" 




> 








T~ x ~ 

X 














< 12 > 


i *-> ■ _ 






i 














- sua 


ween 


epao 


lug J> 








T 





























.02 



.04 



.06 .08 
6, radian 



.10 



.18 



.14 



Figure 9.- Computed and experimental values of variation of angle of twist 
with torque for box beams. Experimental data from reference 8. 
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